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Abstract 

We introduce the strong CP-violation in the framework of AdS/QCD model and calculate the 
electric dipole moments of nucleons as well as the CP-violating pion-nucleon coupling. Our holo- 
graphic estimate of the electric dipole moments gives for the neutron dn = 1.08 x 10^^^ 9e ■ cm, 
which is comparable with previous estimates. We also predict that the electric dipole moment of 
the proton should be precisely the minus of the neutron electric dipole moment, thus leading to a 
new sum rule on the electric dipole moments of baryons. 
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I. INTRODUCTION 



The strong interactions of elementary particles are known to be highly symmetric. The 
most stringently tested global symmetry of strong interactions in the framework of relativis- 
tic quantum field theory is CP, charge conjugation (C) times parity (P), or the time reversal 
symmetry, T. The experimental upper bound on the CP violation comes from the absence 
of the electric dipole moment (EDM) of neutron, 

< 2.9 X 10"^^e-cm (1) 

with 90 % confidence level hi . More precise measurements of the neutron EDM are now 

[1 

under way to improve the current limit by a factor of 50 to 100 [2|. On the other hand, the 
standard theory of the strong interactions, QCD, allows a CP violating term, called the 9 
term, 

Q 

-^QCD 9 ^nvpaG^'^^'G"'^'' , (2) 

where is the field strength tensor of gluons. 
The quark mass term in the QCD Lagrangian, 

£„ = -gl(M,),,g^ + h.c. , (3) 

shifts 6* by a chiral rotation, 9 = 9 + ArgDetM^, which is the physical^strong^P- violation 



angle. The 9 term contributes to the neutron electric dipole moment 



y, y, fl y, I9 



dn = cx 10-^^ 9e-cm, (4) 

where c is a constant of order one, as can be seen from the naive dimensional analysis, 

1 \ / m 



dn-( — ]C^ (5) 

with mjv and niq being the masses of nucleons and quarks, respectively. The 9 parameter 
therefore has to be extremely fine-tuned to be consistent with the experimental data, 

^ + ArgdetM < 10"^ (6) 

Such fine-tuning is known as the strong CP problem and several solutions are proposed to 
solve the strong CP problem |lo| |. 
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Since not only QCD but the electroweak interactions ll| and also the physics beyond 
the standard model contribute to the neutron electric dipole moment, it is quite important 
to estimate the QCD contribution accurately. In this letter we estimate the electric dipole 
moment of nucleons as well as the CP- violating pion-nucleon coupling in holographic models 
of QCD, which have been quite successful in describing the properties of hadrons. We find 
that our holographic estimate of nucleon EDM is comparable with previous results, based on 
lattice calculations Q], current algebra [s], 4,0], chiral perturbation theory or QCD sum 
rule 9|. We also get an interesting sum rule for the nucleon EDM's; the EDM of neutron is 

Q 

opposite to that of proton, dn + dp = 0, which is consistent with the recent lattice result [7|] 
and the estimate from the Light-Front formalism {l2 |. 

The new sum rule on EDM is insensitive to any higher order corrections in l/N^ and 
is also a model-independent prediction of holographic QCD, where baryons are realized as 
instanton solitons. As was shown in Q, the Pauli term in the 5D action should not have any 
f/(l) coupling, since the instanton has only SU(2) nonabelian tails, and thus the anomalous 
magnetic moments of baryons should add up to zero for each flavor multiplets. The same 
should hold for the electric dipole moments, because they are related to the anomalous 
magnetic moments by a U{1) axial rotation. 



II. THE MODEL WITH BARYONS 



A holographic model for spin ^ baryons is constructed for two flavors in 



section, we briefly summarize the model 



14 



In this 



16| , since we will be studying the nucleon EDM 



in the framework of this model, closely following the notations. For the meson sector, we 



take the simplest hard- wall AdS/QCD model as in 17|, [iSj, with the metric 



ds'^ = — {—dz^ + Tj^ydx^dx^^ 



(7) 



where < z < and rj^y = diag(+l, —1, —1, —1). This model captures the important 
aspects of low energy chiral dynamics of light mesons, especially that of pions and vector 
mesons. One should think of this type of holographic models as alternative effective theories 
of strongly coupled field theory like low-energy QCD in the large limit. The theory is 



See Ref. 



15[ for a model of higher spin Regge trajectory. 
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expected to have a classical nature in the large Nc approximation. In conjunction with the 
renormalization group invariance, which is an essential element of quantum field theories, 
the resulting large Nc classical master field should develop a new, dynamically generated, 
space, which corresponds to the energy scale or the renormalization scale of the Wilsonian 
type. The large classical nature and the Wilsonian renormalization group fit together in 
the extra dimension. The low-lying spin ^ baryons such as a proton-neutron isospin doublet 
and its excitations are shown to be naturally realized as 5D Dirac spinors in this picture jl^ . 

Since we are interested in the spin i, isospin | baryons, we introduce (Dirac) spinors 
in our 5 dimensional AdS slice as a holographic realization of spin | baryons.^ However, 
there are two caveats we have to be careful about. The first one is the representation of 
our 5D holographic baryon fields under the chiral symmetry SU{2)l x SU{2)r of QCD, 
which becomes a gauge symmetry in the dual 5D model. We know that the lowest-lying 
4D excitations, the nucleons, form a doublet under the isospin SU{2)j, the diagonal part of 
SU{2)l X SU{2)r, after chiral symmetry breaking, but there is no unique way to assign the 
nucleon charges under the original SU{2)l x SU{2)r symmetry, since the Nambu-Goldstone 



fields can be always multiplied to the nucleon fields |l9l ]. 

For Np = 2 case however, there is an answer. To match the UV anomaly of SU{2)l x 
SU{2)ji X U{1)b from massless chiral quarks {uL^dLj and {uR,dfi) with N^, = 3, there must 
exist massless chiral baryon doublets {pl,^^) and (pr,^^) with the representations (□, 1) 
and (!,□) under SU{2)l x SU [2) p respectively, when the theory is confining but in the 
chirally symmetric false vacuum 20|] . Their baryon charge is Nc = 3 times the quark baryon 
charge. 



Chiral symmetry breaking in the true vacuum introduces a mass coupling for nucleons, 
^xSB ~ -ruN I ] ^ ( PR ur) + h.c. , (8) 



where E = exp(^) G SU{2) is the non-linear group field in the broken phase, which 



rt 



transforms non- linearly asE ^ Ul^U^ under SU{2)l x SU{2)r. As (S) = 1 in the true 
vacuum, the above is invariant under the isospin SU{2)i for which nucleons form a doublet. 
In the symmetric (false) vacuum, E and the isospin singlet the sigma meson (cr) will be 



Note that a 5D spinor has 4-components, like a Dirac spinor in 4D. 
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completed to a linear field X which is bi-fundamental (□, □) under SU{2)l x 5'f/(2)ij, and 
we should have the following term in the theory 



Cm = -9['^ \x(pj,nn)+h.c. , (9) 




to have the coupling in the broken phase {X) ~ AqcdI- 

Based on the above consideration, we find the simplest choice is to introduce two 5D 
spinors Ni and N2, of representation (□, 1) and (1, □) respectively under SU{2)l x SU{2)r 
5D gauge symmetry Upon the Kaluza-Klein (KK) reduction to 4D, the modes from A^^i 
and N2 must include the above-mentioned massless chiral baryon excitations {pi^ni) and 
(PRt'^R) respectively, in the unbroken chiral-symmetric limit. This requirement uniquely 
fixes the IR boundary conditions for A^^i and N2 a^t z = Zm- Then, the natural holographic 
realization of the nucleon mass coupling is to introduce a gauge invariant 5D interaction 

= -gNiXN2 + h.c., (10) 

where X is the bi-fundamental 5D scalar field of (□,□), whose vacuum expectation value 
(VEV) breaks the chiral symmetry in the model. The above coupling will induce the mass 
coupling ([8]) between the would-be massless 4D chiral baryons from A^i and N2 in the chiral- 
symmetry broken vacuum X ~ The coupling strength g must then be fitted to 

reproduce the nucleon mass rriN = 0.94 GeV as the lowest mass eigenvalue. 

The other caveat in our holographic baryon model is the question of chirality in the 5 
dimensional context. Though A^^i {N2) is the holographic dual field to the 4D left(right)- 
handed nucleon operator,^ there is no chirality in 5 dimensions. The 4D chirality is in fact 
encoded in the sign of 5D Dirac mass term 211]. For a positive 5D mass,^ only the right- 
handed component of the 5D spinor survives near the boundary 2; — >■ 0, and this acts as a 
source for the left-handed chiral operator in 4 dimensions. The story is simply reversed for 
the opposite sign case. The magnitude of the 5D mass is given by the AdS/CFT relation 

ml = (A-2f , (11) 



We are using the same notation for both the operator and the state of the 4D chiral baryons without 
confusion. Strictly speaking, iVi is dual to the left-handed chiral operator made of three left-handed 
quarks, and vice versa for N2- 
^ Our convention is £,„ = —m^NN and — {—i)zj^ with j^tpL — 
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where we take A = | for a composite baryon operator of three quarks. Considering the 
anomalous dimension of the composite operator might lead to slightly different results. 

To wrap up the above discussions, our 5 dimensional holographic model for spin ^, isospin 
i baryons in A''^' = 2 sector is given by 



•S'kin = dz dx \/G, 



Sm = j dz j rfx^v^ [-gN^XN2-gN2X^N,] , (12) 

where Dm is the gauge and Lorentz covariant derivative, a/Gs = ^, and the gamma matrices 
in our AdS^ are related to the 4D gamma matrices as = z'-f^ for fi = 0,1,2,3 and 
= —iz'y^. Upon KK reduction to 4D, it is easy to find the eigenmode equations for the 
mass spectrum of 4D spin ^ baryons. Writing Ni{x,z) = fiL{z)BL{x) + fiR{z)Bji{x) and 
similarly for ^"2(3^, z) = f2L{z)BL{x) + f2R{z)Bii{x), where -Bl,_r are the components of the 
4D eigenmode spinor B = [B^, B^)^ with mass mjv to be determined, we have 




) 


(;::) 






)l 


flR \ 








TTlN \ , (13) 



with A (= I in our case) in general. As mentioned before, the existence of the chiral zero 
modes when (X) = requires us the IR boundary condition fiR^Zm) = f2L{,Zm) = 0. In 
the meson sector the best fit was found for (X) = \{ rnqZ + az'^) with niq = 2.34 MeV, 
a = (311 MeV)^ and the IR cut-off z^ = (330 MeV)"^ [3, Q. Then, the only remaining 
parameter of the theory is the dimensionless coupling g, which was found to he g = 9.18 
in |l^ to reproduce = 0.94 GeV as a lowest mass eigenvalue. 

In [1^, which discussed nucleons in the top-down Sakai-Sugimoto model, an important 
new operator, responsible for anomalous magnetic dipole moments, among others, was iden- 
tified, which will also be important for our analysis of CP-violating electric dipole moments. 
Experimentally, the proton magnetic moment is /ip = 2.8fiN = I^n + l-SyUAr where the Nu- 
clear Magneton /iat = is from the minimal coupling of charge 1, while the neutron has 

= — l.S/iAT. It is clear that the anomalous piece has a structure of isospin doublet without 
U{1)b charge, as shown in |l3]. Since electromagnetic coupling is a sum of the isospin and 
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U{1)b, Q = I-B + /3, the operator responsible for the anomalous magnetic moment does not 
include U{1)b- 

It is easy to find the corresponding 5D operator which induces the desired anomalous 
magnetic moment upon 4D reduction, 

C^ipoie = iD [N^T''''{Fl)mnN^ - N2T^'''{Fr)mnN2] , (14) 

where Fl^r are the field strengths of Al^r and D is a real parameter that must be fixed 
to reproduce //anomalous = 1.8/iAr.^ Note that Al^r in our model do not include U{1)b, as 
required. In fact, this absence of U{1)b in the 5D effective operator has its origin from 
the fact that baryons in AdS/QCD arise as instantonic solitons of small size in our (4+1)- 
dim gauge theory of A^^r. As we treat them as point-like with the effective fields A^^i and 
iV2, we need to take into account their long-range instanton tail of Al^r, and the operator 
exactly sources those tail profiles. In a more complete description, the coefficient D 
would be determined by the stabilized size of small instanton-solitons, whereas here it is a 
fitted parameter against experiments. Since instanton profiles are purely non-abelian, the 
operator (IT^ should not include U{1)bi which leads to a model-independent sum rule for 
any quantities derived from the operator such as the anomalous magnetic moments. 

III. PHYSICS OF VACUUM ALIGNMENT 

The strong CP-violation in QCD can be introduced either in the vacuum ^-angle or as 
complex phases in the quark mass matrix Mq defined as 

= -gi(M,)i,g]j + h.c. . (15) 

As is well-known, the anomalous axial U{1)a rotation can shift the 5-angle to zero, making 
the strong CP-violation appear only in Mq or vice versa. The physical strong CP-violation 
angle is^ 

^ = ^ + ArgDet(Mg) . (16) 

^ In |l3|. D was reliably predicted from string theory, which explains /^anomalous = \.%^n quite well. But in 
our bottom-up model it is a fitting parameter, though the relative minus sign in (I14p is dictated by 4D 
parity invariance, as easily seen from the fact that 5D Dirac mass term flips its sign under 4D parity. 

^ The 0-angle is normalized as Cg — 3^Tr(i^ A F) with J Tt{F A F) — 1 ior a, single QCD instanton. 
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Presumably, the ^-angle can be described in the holographic model by a kind of axion |22| . 
but for our purpose it is much more convenient to work in the frame where the strong CP- 
violation appears in Mq only. This is because in our AdS/QCD model, Mq is easily described 
by the non-normalizable mode of the scalar X, 

{X{z)) = \MqZ + h:z\ (17) 

where S is the quark bi-linear condensate (S)*-' = — 

In the CP-symmetric case of Mq = diag(m„, m^^) with small real masses m^^ <^ Aqcd, 
it looks natural to have the bi-quark condensate proportional to identity matrix^ 

S ~ ^Qcola , (18) 

as was used before and will be shown rigorously in a moment. However, as we turn on small 
CP-violation as phases ai^2 of quark mass Mq = diag(e*"^m„, e^^^m^), it is not obvious which 
polarization S will take in the SU (2) space of its moduli, S = A'^U with U G SU (2) and A 
is a scale proportional to Aqcd, the intrinsic scale of QCD. This is the problem of vacuum 
alignment, first considered by Dashen [23I ]. 

The vacuum moduli of S form a coset space SU{2)l x SU{2)r/ SU{2)i. Since the quark 
mass Mq breaks the chiral symmetry, it lifts the moduli space by inducing a potential, 

F™ = -Tr(M,Et)+h.c. , (19) 

which can be inferred from ( fT5l) . With our diagonal mass matrix M^, it is natural to make 
(S) align along its diagonal direction by a SU{2)l x SU{2)^ chiral transformation, that is, 

(gix g \ 
(20) 
e-^y 

to get the vacuum potential Vm = — 2m„ A^ cos(x — ai) — 2mdA^ cos(x -|- 0:2). For small 
ai^2, related to ~ 10^^, as will be seen in a moment, we expand the above cosines to get 
Vm/A^ ^ {rriu + md)x^ — 2(m„ai — mda2)x + const. The potential minimizes at 

rriu + rrid 



Isospin violating effects from niu 7^ ind to the bi-quark condensate can be neglected to first order in m^, 



as we are interested in the results that are first order in rriq and 0. 
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as the aligned vacuum polarization in the presence of the strong CP-violation. It also 
confirms that identity matrix is the ground state for the CP-symmetric case, a\^2 = . 

For our analysis, it is more convenient to take a non-anomalous SU{2)a axial rotation 
to remove the phase angle x in S.^ This brings the quark mass angles to ai —>■ ai — 



i^^i±^ and a2 ^ a2 + x = i^^i±^2)lll^. 
9 = («! + 02)) "we finally end up with the quark mass matrix 



X = '."1 a2)md Q, ^ I 2; — (ai Q2jm„ ^ gjnce the physical CP- violation angle is 



M,= \ " ^g,^^. I and (S) = al . (22) 

m^e 

Therefore, in our AdS/QCD model the strong CP- violation is easily encoded in the VEV of 
X as {X{z)) = I (MgZ + (S) z^) with Mg and (S) given in ([22]). This is the starting point 
for our holographic analysis of the strong CP-violation. 

Note that we introduce strong CP- violation only through the VEV of X, while keeping the 
AdS/QCD theory itself, namely the various couplings and all other parameters of AdS/QCD, 
CP-conserving, because the QCD dynamics is CP-conserving when the angle 6 vanishes, and 
the angle can be absorbed into "external" quark mass parameters. 

As we are interested in the effects of the strong CP violation in the first order of 6 (and 
mu,d), we expand the phase exponents in Mg to get 

(m \ 
\ +i9ml = + i5Mg , (23) 
m,, J 

where = is the reduced mass. The CP-violating axial mass i SMg is propor- 

tional to the identity matrix ^, which could also be argued from the vacuum stability 3|, |4|. 
This in fact agrees with the previously known chiral Lagrangian with the strong CP- violation. 
For simplicity, we will take an isospin symmetric quark mass m„ = rrid = m = 2m, then the 
whole {X) is proportional to the identity matrix, which we write as 



1 i - ' 

-imz + az^) H — mO z 
2^ ^ 4 



v{z)l = [vq{z) + i5v{z)]l , (24) 



where Vq{z) is the dominant CP-conserving piece while small i5v{z) = i^z breaks the CP 
symmetry. 



^ Note that this global rotation also rotates Mq simultaneously. When we discussed the lifted moduli of 
SU{2) above, we fixed Mg. In other words, only the relative SU{2)a angle between Mq and E is physical 
and we can always go to the frame where E ~ 1, as the coset space is homogeneous. 

^ The angle 9 should be proportional to 75, which is suppressed here, since we are using the chiral basis. 
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IV. THE NEUTRON ELECTRIC DIPOLE MOMENT 



The effects of the strong CP-violation in the VEV of X are generically of two kinds in the 
holographic models; it modifies the wave functions of 4D excitations along the fifth direction 
if the corresponding 5D field couples to (X), or it affects various couplings of 4D interactions 
after integrating over the fifth direction if the 5D coupling involves X. In consideration of the 
latter, we should also take into account the former as well. The CP-violating pion-nucleon 
coupling g-^NN that we will calculate in section |V] is such an example. 

The operator f|T^ will be responsible for the electric dipole moment of the nucleons, as 
well as the magnetic dipole moment. This is because an axial U{1)a rotation of the magnetic 
dipole moment is precisely the electric dipole moment, and the strong CP- violation can be 
described by the non-zero U{1)a phases of the quark mass matrix, so that its effect will 
appear as a small U{1)a rotation of the magnetic dipole moment. We stress that the 5D 
operator f|T^ itself is CP-conserving, and we don't include explicit CP- violating 5D operator 
in the model. After introducing the CP-violation in (X) in the otherwise CP-conserving 
theory, additional CP-violating 5D operators may be generated by 5D loop effects. However, 
in the spirit of AdS/QCD we keep the tree-level analysis only, systematically ignoring 
corrections. Therefore, we expect our analysis captures a leading Nc contribution only 
except the sum rule, dp + = 0, whose origin goes beyond the l/N^. approximation. The 
CP-violation introduced by (X) affects the 5D wave functions /(i,2)(L.i?,) of the (lowest) 
nucleon state through the mass coupling 

/:m = -9Ni{X)N2 + h.c. , (25) 

and the resulting 4D operator of nucleons after integrating over the fifth coordinate with 
these wave functions will have a small CP-violating counterpart of the magnetic dipole 
moment, which is the electric dipole moment we are interested in. 

To be explicit, the CP-violating (X) = v{z)l enters the equation for the nucleon wave 
function Xi = fiL{z)BL{x) + fiR{z)BR{x) and X2 = f2L{z)BL{x) + f2R{z)BR{x), where 
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Bl,r{x) is the 4D nucleoli field, as follows 





1 


^ flL 












^ flR 




I 


[ f2R 



-niN 




(26) 



where A = | is the scaling dimension of the 4D baryon operator and mjv = 0.94 GeV is 
the nucleon mass. To get the standard 4D kinetic term for B{x), we have to normalize the 
eigen-functions 



dz 



1 



-A\hL\'+\f2Ln 



dz 



(27) 



We are interested in the effects of small CP-violating imaginary part in v{z) = vo{z) +iSv{z) 
to linear order. The CP-conserving zeroth order eigenfunctions /(i 2)(l r) with vo{z), obtained 
numerically in jl^, are real functions with an important property 



/: 



(0) 
IL 



AO) 

JIR 




(28) 



which can also be checked explicitly from ( l26i) . This is a simple consequence of 4D parity. 
The precise parity transformation in our model that leaves the 5D action invariant is 



[Xq,X,Z) < — > [Xo,-X,Z) 



(aIAl^ai) 




(29) 



We see that the (lowest) nucleon state fl28|) is parity even as expected. Using this and the 
explicit form of 6v{z) = mO z/A, it is not difficult to find that the first order effect of i6v{z) 
on the wave functions is a simple phase factor in the eigenfunctions /(i,2){l,_r), 




AO) 

JlL 

AO) 

J2L 




JI3 



AO) 

JlR 

AO) 

J2R 



with 



{a -(3) 



9 f 5v{z) \ ^ g m 
rriN \ z J ArriN 



(30) 



(31) 
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without affecting the mass mjv = 0.94 GeV. Note that the common phase on fi and Jr is not 
physical as it wouldn't appear in any bilinear operators Nr'^^'^"'N, while the above difference 
in phases f l31l) is a physical CP- violating effect. It is also possible to work in the frame where 
all wave functions are real, and the mass uin instead gets a phase mj\[ e*''"~^^m7v/° 
although we will stick to the previous description for a while. We will get back to the frame 
of complex mass and real wave functions later, when we argue the generality of our result 
against other possible higher dimensional operators. 

Inserting the expansion Ni = fiL{,z)BL{x) + fiR{z)BR{x) and N2 = f2L{.z)BL{x) + 
f2R{z)BR{x) into our relevant 5D operator (JT4l) . 



S-dipoie = iD j d^x j dz^/Gl [N^V^''{Fl)mnN^ - N^V'''' {Fr)mnN2] , (32) 

it is straightforward to obtain the resulting 4D dipole moment operators. To read off the 
electromagnetic coupling, recall that Q = h + ■ ■ ■ with the isospin Is = (tg + 1^) , so that we 
can simply replace Al^r by = eA^^t^ and Ar = eA'^'^ts with = ^cts to extract couplings 
to the electromagnetic vector potential A*^™. According to the AdS/CFT correspondence, 
the external vector potential is simply a non-normalizable mode of the corresponding 5D 
gauge field, which turns out to be a simple constant mode over the fifth direction. 
A quick calculation results in the following 4D dipole moment operators, 



-^^magnetic - ( T^Ba^^a^B j F^™ , 



Electric = 4 (^-^^^^'^7^35^ F^^ , (33) 

where acts on the isospin index of the nucleon doublet B = {p, n)^, and the anomalous 
magnetic dipole moment fi'^" as well as our CP-violating electric dipole moment de are given 
by 

" 1 



/zr = e-(-l)-(2Z}) / dz 

Jo 

de = e-{a-(3)- {2D) / dz 



_ AO) AO) 







}_ AO) AO) 

^sJlL J2L 



(34) 



In the general complex mass case, the mass rriAr in the first equation in ((26|) should be replaced by m 



N- 



Our convention is ct^'' = ij^^" = §[7'^, 7''] with 7^ = ( ) and 7^ = [ ^ | • Here cr^ = (1, -a) 



\ -1 



and cr^ = (1, a). 
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where factor in the integrals is traced back to ^/Gl = and the curved space gamma 
matrix carries an extra factor of z since the inverse vielbein = . The coefficient D 
must be determined by comparing with the experiments jJ^" = l.Sfijy = but for our 

purpose of obtaining de, we have a model- independent prediction for the ratio 

-(«-/?) = -!— ^ , (35) 



between the electric dipole moment and the anomalous magnetic dipole moment of nucleons. 
This gives us 

4 = -1.8,„.f^»=-i|*(^)(^)». (36) 

ArriN o yrriN J \mN J 

Using niN = 0.94 GeV and GeV"^ = 0.197 x 10~^^ cm, we have 

4 = -0.25 g X [^^) X lO-i*^ e (e ■ cm) , (37) 

in units of (e ■ cm). This is the main result of our paper. 

For a model in (l7| with Zm = (330MeV)"\ a = (311 MeV)^ and m = 2.34 MeV, we 
have g = 9.18 to have the correct nucleon mass eigenvalue m^v = 0.94 GeV [l4], and this 
predicts that 

4 = -1.08 X 10"^^ 9 (e ■ cm) . (38) 

This means that the neutron electric dipole moment is 4 = —4 = +1.08 x 10"^^ 6 (e-cm). 

We end this section by pointing out a new feature that our result predicts for the first 
time. The isospin structure in the operator ( l33i) without U{1)b component tells us that the 
proton electric dipole moment dp should be equal in magnitude but minus to the neutron 
electric dipole moment, that is, dp = —dn = d^.. It would be very interesting to test this 
prediction in the future experiments, which would confirm the validity of our 5D model of 
nucleons. 

V. CP- VIOLATING PION-NUCLEON COUPLING 

We calculate the CP-violating component in the coupling of the pions to the nucleons 

'Ccpodd = g-KNN {Ba B) -Ti , (39) 

where B = {p, n)'^ is the nucleon iso-doublet and a acts on the isospin index, while the usual 
CP-even interaction is defined as 

^CPeven = WitNN { B a 'J^ B ) ■ TV . (40) 
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As we have the nucleon wave functions /(i,2)(l,_r) in the previous section, what remains is 
to identify the pion profile along the fifth direction, so that the above 4D interactions can 
be read off by integrating out 5D operators over the fifth dimension. We will pay a special 
attention to how the strong CP-violation we introduce in {X) = [vq{z) + i6v{z)]l induces 
a mixing between the CP-even pions and the isoscalar mesons from X, since this will also 
contribute to (139|) eventually. 

Normally, the pions are identified as the axial SU{2)a phase angles of the chiral conden- 
sate, that is, X ~ (X) e*^(^)^(^) with a 5D profile f{z)}^ In the AdS/QCD set-up, since we 
elevate the axial SU (2) a to a gauge symmetry, this is no longer a gauge invariant statement, 
and we should consider the broken 5D SU{2)a gauge field simultaneously with the phase 
angles of (X) which could be viewed as eaten Goldstone bosons by the (5D) Higgs mecha- 
nism. In our analysis, we choose to work in the unitary gauge, or equivalently a 5D version 
of the -R^-gauge in ^ — » oo limit, which will become clear in a moment. The advantage 
of this gauge-fixing is an unambiguous diagonalization of the kinetic terms of the relevant 
fields, which makes finding equations of motion a lot easier than other gauge choices. 

A modification appearing with the strong CP-violation turns out to be that not only the 
phases but also the modulus of the chiral condensate takes part in the pion wave-functions; 
in other words, expanding X ~ {X) e"'^^*^ with P and Q being Hermitian, there is a small 
CP-violating mixing between the modulus Q and the SU{2)a gauge field in addition to 
the usual mixing of P and Az}^ {A^ is the component of the SU{2)a gauge field along the 
fifth direction z). This implies that the resulting pion wave- function along the fifth direction 
involves a small Q-component as well as the dominant Az and P part. The CP- violating 
Q-component in the pion wave-function will then contribute to the CP-violating coupling 
to the nucleons in fl39|) . Of course, the operator fl39l) also receives contributions from the 
CP- violating phase (a — P) in the nucleon wave-functions /(i,2)(L,i?,) that we identify in the 
previous section. 

This corresponds only to the phase of the chiral condensate, without affecting the current quark mass. 
The reason is that f{z) is normalizable and decays to zero in UV region, and the IR profile of X encodes 
the chiral condensate, while its non-normalizable UV mode contains the current quark mass. Therefore, 
the normalizable pions with /(2;)-profile do not contain the phase of the current quark mass. 
^•^ There does not appear any CP-violating mixing between P and the vector SU{2)v since, in our set-up, 
{X) is proportional to identity matrix (X) — v{z)l — [vq{z) +idv{z)]l, and the SU{2)v remains unbroken. 
This may not hold in general for an isospin asymmetric quark mass niu ^ ra^. 
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5^5] = / d''xdz^/G^TT{ - 



It is rather straight-forward but tedious to perform the above mentioned procedure to 
obtain the 5D profile of the pions. Defining the axial vector field hj A = ^{A^ ~ Aji) and 
expanding X = {X) e"'3+*^ = [t>o + i6v]e~'^~^''^ up to the quadratic order, we have 

I - ^(Fa)m^(Fa)^^^ + {vof{2AM - OmP) {2A'' - d^'p) 

- 4(^{d,Vo)Sv-Vo{dJv)^Q{2A^-d'P^^ , (41) 

where Fa is the field strength of A and the last line is the CP-violating mixing in the linear 
order of 6v. In the above, we should remember that \/G^ = and = g^^d^ = —z^d^, 
A^ = —z^Az- As we are interested in the effects from the small CP- violating bv in the last 
line, we first solve the top line for A^ and P to find the zeroth order profile of the pion, and 
then treat the last line as a perturbation to calculate a small induced Q-component in the 
pion wave-function. 

The kinetic terms are mixing among A^, A^ and P in the first line of (1411) . and to remove 
those mixing terms, it is convenient to introduce a 5D version of gauge- fixing term jl^ . 

n 2 

(42) 



Sa f = / Ct^xdz - 

2i J z 



where everything is written in the flat metric basis, ds^ = rj^j^dx^dx^ — dz'^. We then have 
separation between and [A^, P), and, since the pions arise from the latter, we focus on 
the [A^, P) sector only from now on. We further take —>■ oo limit to work in the unitary 
gauge for simplicity. This means that we need to impose a constraint for the propagating 
fields in such a limit, 

SO that P can be replaced by Az in the action. We then solve the resulting action for A^ to 
find the pion profile. The resulting relevant action is 



+ ^^{{d,Qf - {dzQ?) --A{dzv,f - ^)g2 - \v,{dzV,)QdzQ 

^3 



z 



+ l((a.„M,,-,,o(ai.))Q(2A.-4(^a(^)))l . (44) 
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where we first solve the top hne and then solve for Q from the last two lines. 

The equation of motion can be easily derived from the above. As we are looking for the 
pion profile, we put = it{x)A{z) and Q = 7r(x) q{z) with q{z) being expected to be small 
and linear in 5v oi 9 . The 4D pion field tx{x) satisfies 9^9^7r(x) = — m^7r(a;) where m^r is 
the pion mass that must be determined from the eigenvalue equation we will specify in a 
moment. The equation for A{z) turns out to have a nice factorization structure, due to the 
underlying gauge invariance. To be explicit, defining 

B^2A-d^ {^^^—,d, {-) ) , (45) 



-2^|(t;o)- \z 

we find the equation for A{z) becomes a equation for B{z) only. 

.3 .,..,2 .\ 4^2(^^)2 




Z2 



B = -mlB , (46) 



with the IR boundary condition B{zm) = 0. The normalizability at UV, z — > 0, determines 
the pion mass rriT^ as the eigenvalue of the above equation. For example, the parameters 



in [17|, Zni = (330MeV)-\ = Att^, and vo{z) = \{mz + oz^) with m = 2.34 MeV, a = 
(311 MeV)^, indeed give us m^r = 140 MeV as the lowest eigenvalue The normalization 
of B{z) should be determined such that the 4D action of it{x) after integrating over the 



z-direction must have the form 



d xTt 



{d,n){d^n)-min' . (47) 



Looking at the last term in the first line of (jS]), which produces the mass term 







""rf^M!^^. / d^xTrin^] , (48) 



we see that the normalization of B{z) must be fixed by 



Once we solve B{z), we then find A{z) via solving (HSj) with the IR boundary condition 
A{zm) = and the normalizability in the UV region, z — > 0. Finally, p{z) in P = p{z)7i{x) 
is determined through the constraint fH3l) . 



One can also easily check that the Gell-Mann-Oakes-Renner relation holds to a good degree in our model 
by calculating the pion mass for different values of the quark mass. 
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Having obtained the zeroth order profile A{z) (and p{z)) for tlie pion, we next solve q{z) 
in Q = q{z)TT{x) that is induced by the CP- violating mixing in the last line of dS]). Note 
that the effect of this mixing to A{z) and p{z) is of second order in Sv, and can be neglected. 
The equation for q{z) is written as 



z" 



z 



3 ^ id.vof - 



(51) 



where the second line is the source for q{z) induced by 6v{z) 



md 



z. We should impose the 



IR boundary condition qi^z^) = and the UV normalizability as well. 

After getting the pion profile = A{z)7r{x), P = p{z)7r{x) and Q = q{z)n{x), we see 
that it is straightforward to insert them together with the nucleon wave-functions /(i,2){l,_r) 
into our 5D action of holographic baryons, and to read off the couplings between pions and 
the nucleons. Writing the resulting 4D interaction terms as 



(2) 



TT 



(52) 



where we expand vr = vf ■ | , the last term is the CP- violating pion- nucleon coupling we are 
interested in, with the coupling strength q-^nn given by 



g-KNN 



-{a — (3) dz 
Jo 



iA{z)fS!{z)fi^iz) + ^M^)p{z) [{fS!{z)r + ifS!i^)f 



1 



Vo{z)q{z) + 6v{z)p{z) ) {fiS{z)r - if^'^'iz)) 



ciO), 



(53) 



where {a — P) = ^j^O and 5v{z) = ^z as given before. This is the primary result of this 
section. Since q{z) is proportional to 9 via 6v, our g^NN is proportional to the strong CP 
angle 9 as expected. 

For the model of Zm = (330 MeV)"^ and vo{z) = \{mz + gz^) with m = 2.34 MeV, 
o = (311 MeV)^, we have 5^ = 9.18 and our numerical result for g-^NN is 



gnNN 



+0.017^ , 



(54) 



which is about half from the previous estimate in 



17 



For completeness, we also give the expressions for the CP-conserving coupling constants 



y-KNN 



(2) 
9-kNN 



dz 



-,A{z)fi'^{z)f^^>{z) + ^Vo{z)p{z) {fl^>{z)r + {f^^{z)f 



p(0). 



p(0)/ 



-(2m7v) -D dz 



A{z) ifS!iz)r + if^'^Hz)) 



(0), 



(55) 



Note that the coupling term with gf^^j^ may turn into the (yf^^^- coupling by using the 
equation of motion for B, and the total pion-nucleon coupling is q^nn = O^nn + Q^nn- 
Recall also that D should be fitted against the anomalous magnetic dipole moments of the 
proton and neutron by ([3] 



dz 



±_ AO) AO) 



u 



1.8 



2m N 



Numerical values for the above model are (y'^jvAf 

cxp 



-3.64, g^;l^ 



(56) 
-25.6, 



—22.0 with gT,NN 

which is about twice of the experiment g^^^^ = 13.6. There is however much chance to 
improve several crude simplifications in constructing our model, such as the simple AdS- 
slicing, neglecting anomalous dimension of the three-quark baryon operators, the possibility 
of a different IR-cutofF, etc. 24i\ . 



VI. DISCUSSION AND CONCLUSION 



We have studied the CP-violating effects in strong interactions, induced by the QCD 
6'-term, in a bottom-up model of holographic baryons [14]. We find that the model of 
holographic baryons is quite useful in studying the CP- violating effects in the baryon sector 
such as the electric dipole moments of nucleon or the CP-violating couplings of pions to 
nucleons, since the CP-violating effects are easily incorporated in the bulk scalar field X, 
whose coupling to bulk spinors naturally leads to the calculable CP-violating effects in the 
baryon sector. 

A Pauli term is added to the holographic model of baryons [l^ to estimate the electric 
dipole moments of nucleons. We have fixed the coefficient of the Pauli term by fitting the 
anomalous magnetic moments of nucleons, though in top-down holographic models such as 
Sakai-Sugimoto model [25| the coefficient can be calculated reliably due to the fact that the 



baryons are realized as instanton solitons in holographic QCD 
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26|. However, we stress 
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that our result for the electric dipole moment of nucleons is model-independent, since it 
relies only on the ratio between the electric dipole moment and the anomalous magnetic 
moment of nucleons, d^/ > which is independent of the coefficient of the 5D Pauli term, 
as shown in fl5^ . 

In fact, we can argue further that our result is universal against any other possible higher 
dimensional operators in 5D. Note that in solving (126!) . we can alternatively go to the 
frame where the nucleon wave-functions /(i,2)(L,i?) are real, whereas the strong CP-violation 
is encoded instead in the phase of the nucleon mass m^v — > e'^^°'~^''mj\j. Upon reduction 
to 4D, all the resulting 4D operators obtained from real wave functions /(i,2)(L,i?), including 
anomalous magnetic dipole moment operator, will be CP-conserving except the nucleon mass 
term with the complex axial mass. Performing an axial U{1)a rotation back to make the 
nucleon mass real, this will induce the electric dipole moment operator from the anomalous 
magnetic moment operator as in (133!) with the universal prediction for the ratio 

.(a- f3) = --—e, (57) 



without any regard to the details of the 5D model. 

Though our model is bottom-up, the bulk spinors in our model should be interpreted 
as effective fields of the instanton solitons as in the top-down model, because the mesonic 
sector of our model, being a 5D gauge theory with a Chern-Simons term, naturally supports 



instanton solitons, whose Wilson line can be identified as skyrmions 
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291 ]. As shown 



m yjij, the Pauli term therefore should not contain the U{1) gauge fields to reproduce the 
correct behavior of non-abelian instanton solitons at long distances, which leads to important 
sum rules for both the electric dipole moments and the anomalous magnetic moments of 
nucleons, which are model- independent and not subject to any l/N^ corrections, 

dp + dn = 0, + /ir = . (58) 

Combining our holographic estimate of the neutron electric dipole moment, given as 
dn = 1.08 X 10^^^ 6e ■ cm, with the current experimental bound on the neutron EDM [l|, 
we obtain a limit on the 6 angle reliable up to 30% coming from 1/Nc corrections, 

1^1 < 3 X 10"^°, (59) 

which is somewhat stronger than the previous bounds. 
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